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A major research and technology program in Probabilistic Structural Analysis 
Methods (PSAM) is currently being sponsored by the NASA Lewis Research Center 
with Southwest Research Institute as the prime contractor [1,2]. This program 
is motivated by the need to accurately predict structural response in an environment 
where the loadings, the material properties, and even the structure may be 
considered random. The heart of PSAM is a software package which combines 
advanced structural analysis codes with a fast probability integration (FPI) 
algorithm for the efficient calculation of stochastic structural response. 


The most common structural analysis techniques in use today are finite 
element and boundary element methods. These techniques permit highly accurate 
three-dimensional modeling of structural geometry, thermal and mechanical load 
environment, and material properties. Unfortunately, these techniques can also 
be expensive and time-consuming. In the early stages of design, when geometries, 
loads, and materials are only tentative, it is not practical to assemble an 
exact, comprehensive three-dimensional model for every critical component. What 
Is needed, typically, is some means of estimating the nature and rough magnitude 
of stresses, displacements, natural frequencies, etc. Of further value is the 
identification of those design parameters which exert the most influence on the 
total system performance, so that further design evolution is more efficient. 

These needs are addressed in the PSAM software by the Probabilistic 
Approximate Analysis Methods (PAAM) module. The basic Idea of PAAM is simple: 
make an approximate calculation of system response, Including calculation of the 
associated probabilities, with minimal computation time and cost, based on a 
simplified representation of the geometry, loads, and material. The deterministic 
solution resulting should give a reasonable and realistic description of 
performance - limiting system responses, although some error will be inevitable. 
If the simple model has correctly captured the basic mechanics of the system, 
however, including the proper functional dependence of stress, frequency, etc. 
on design parameters, then the response sensitivities calculated may be of 
significantly higher accuracy. In other words, the calculated probabilistic 
distribution of the response variable may be in significant error only by some 
offset of the mean value. 


Three factors make up the ’’approximate" analysis approach. The first Is a 
simplified representation of the part geometry. Complex three-dimensional shapes 
are replaced by simple beams, shells, etc., with relatively few descriptive 
parameters. The second factor is a similar simplification of the applied 
mechanical and thermal loads. Point loads, uniform fields, and linear or parabolic 
distributions are used to describe these quantities . Often a complex load 
environment can be approximated through linear superposition of the simple 
descriptors. The third factor In the PAAM methodology is the solution technique 
itself. Complex numerical techniques with many degrees of freedom are replaced 
by simpler computational schemes. In some cases, it is possible to construct 
a simple mechanics - of -materials model, often in two dimensions. In other cases, 
more sophisticated closed-form solutions can be derived or (preferably) adapted 
from previous research results in the literature. These maybe based on elasticity 
or may employ approximate energy methods . 
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The PAAM code employs the fast probability integration (FPI) algorithm of 
Wu and Wirsching [3], which has been demonstrated to consistently provide fast 
and accurate estimates of point probabilities for typical engineering response 
functions. This algorithm identifies a most-probable -point (or design point) 
on the response surface, establishes a quadratic polynomial approximation to the 
response function at that point, and then transforms the quadratic form into a 
linear one. Input variables are defined using probability distributions. An 
optimization routine is employed to approximate non-normal variates as equivalent 
three-parameter normals, thereby approximating the limit state as linear in 
normally distributed design factors. Solution of this simplified problem is 
straightforward. The FPI is repeatedly applied to a number of limit states to 
compute the cumulative distribution function (CDF) , 

The PSAM code is being demonstrated and validated by analyzing four 
representative critical components in the current Space Shuttle Main Engine 
(SSME) . The same four components are considered in the PAAM code, which by its 
very nature Is component- specific . It should be noted, however, that the 
architecture of the PAAM software permits new closed-form expressions for these 
or other components to be Installed and evaluated quickly. In this presentation 
one of these demonstration components, the liquid oxygen (LOX) post, Is used to 
Illustrate typical PAAM solution strategies and results. 

Two different simple models are used to analyze the LOX post. The steady 
temperature and pressure gradients between the Internal and external radii are 
evaluated by modeling the post as an ax i symmetric thick cylinder. This is a 
classic elasticity problem, and exact solutions are readily available in standard 
textbooks. Other loadings, including the transverse fluid flow, are considered 
by modeling the post as a beam with hollow circular cross-section and elastically 
restrained ends (i.e., a mechanics-of -materials approach). Note how this model 
is used to solve a vibration problem. The basic solutions for free vibration 
of fixed- fixed and fixed-pinned uniform beams are available in handbooks. 
Solutions for a limited number of intermediate end conditions are available in 
tabular form. These results can be fitted with a simple empirical relationship 
to estimate the frequency factor for any mode and any end condition. Subsequent 
calculations of mode shapes and stresses for forced vibration follow directly 
from simple equations. 

Several sample PAAM/FPI calculations have been performed in order to validate 
and further demonstrate the approximate methods approach. Results are given in 
terms of cumulative distribution functions for critical response variables. FPI 
and Monte Carlo calculations are shown to agree closely. Also determined are 
the sensitivity factors of the random variables. These sensitivity factors are 
a function not only of the mathematical form of the equation, but also of the 
distribution types and coefficients of variation of each random variable. They 
are of particular value for guiding improvements in the analytical scheme or 
component design. 
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AN OVERVIEW OF PROBABILISTIC STRUCTURAL ANALYSIS METHODS 


0 PSAM LINKS STRUCTURAL ANALYSIS 
WITH FAST PROBABILITY INTEGRATION 
FOR EFFICIENT CALCULATION OF 
STOCHASTIC STRUCTURAL RESPONSE 
WHEN LOADS, MATERIALS, AND 
GEOMETRY ARE RANDOM 

0 FEM AND BEM ARE MOST COMMON 
STRUCTURAL ANALYSIS TECHNIQUES 

- DETAILED MODELS 

- ACCURATE 3-D SOLUTIONS 

- EXPENSIVE, TIME-CONSUMING 

- NOT APPROPRIATE FOR 
EARLY DESIGN/ANALYSIS 

o Need a rapid, inexpensive estimate 

OF STRUCTURAL RESPONSE WHICH 
INCLUDES PROBABILISTIC INFORMATION 



0 "PROBABILISTIC APPROXIMATE ANALYSIS METHODS" (PAAM) 


WHAT IS AN "APPROXIMATE METHOD" ? 


o Simplified representations of 

- COMPONENT GEOMETRY 

- APPLIED LOADS 

- MATERIAL PROPERTIES 

o Rapid, efficient solution scheme 

- MECHANICS-OF-MATERIALS 

- SIMPLE ENERGY METHODS 

- CLOSED FORM OR TABULATED 

LITERATURE SOLUTIONS 

o Calculated mean value of response 

FUNCTION CONTAINS SOME ERROR 

o If basic mechanics are correct, 

RESPONSE DISTRIBUTIONS AND 
PARAMETER SENSITIVITIES WILL 
BE ACCEPTABLY ACCURATE 





V 
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PAAM CODE USES FAST PROBABILITY INTEGRATION ALGORITHM 


0 Wu-WlRSCHING FPI ALGORITHM 


o Shared by other PSAM modules 
o Identify MPP on response surface 

FOR GIVEN PROBABILITY LEVEL 


jOFHT PROB. OENSrTY IfX) 


0 


0 


0 


Establish quadratic approximation 

TO THE RESPONSE FUNCTION AT MPP 
AND TRANSFORM TO A LINEAR FORM 

Input variables defined using 

PROBABILITY DISTRIBUTIONS 

Approximate non-normal variates 

AS EQUIVALENT 3-PARAMETER NORMALS 



o Solve the resulting simple problem 
(linear with normal design factors) 


PROB. t T « 3001 - / ... J»(X]X (PHOT. OF FAJLURE) 
0 . 


PAAM CODE DEMONSTRATED WITH SSME LOX POST 


o One of four SSME 

demonstration components 

0 A SLENDER, hollow tube 
FIXED AT BOTH ENDS 

o Important loadings include 

- AT AND Ap ACROSS TUBE WALL 

- STEADY AND DYNAMIC LOADS DUE 

TO TRANSVERSE FLUID FLOW 

0 TWO DIFFERENT PAAM MODELS SHOWN 

- THICK CYLINDER MODEL 

- BEAM MODEL 

o Although these equations 

ARE COMPONENT-SPECIFIC, 

the basic PAAM algorithm 

IS COMPONENT-INDEPENDENT 


ROW 



ROW 1) 
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THICK CYLINDER MODEL USES KNOWN ELASTICITY SOLUTIONS 


o Thick cylinder subjected to 

DIFFERENT INTERNAL AND EXTERNAL 
PRESSURES AND TEMPERATURES 

o Exact closed-form solutions 

AVAILABLE IN STANDARD TEXTBOOKS 

o Consider the hoop stress: 

- PRESSURE SOLUTION 

(r 0 /r) ? (p„- P.) p,~ p 0 k 2 

°° —\ — + ~i^r 

- THERMAL SOLUTION 

A<jn ?+ I l-1n(r 0 /r) 

Jt*-i 

g/iXT^TJ 
??(l-v) 
k - r 0 /r, 



BEAM MODEL USES A MECHANICS-OF-MATERIALS APPROACH 


o Slender beam with hollow cross-section 

- FIXED AT ONE END 

- ELASTICALLY RESTRAINED AT THE OTHER 

o Consider the vibration problem 
o Basic equations in handbooks 

IN TERMS OF A FREQUENCY FACTOR K 

o Lambda depends on end fixity and mode number 

- CLOSED FORM EXPRESSIONS FOR FIXED, PINNED 

- TABULATED SOLUTIONS FOR ELASTIC RESTRAINT 

- FIT THESE WITH AN EMPIRICAL EXPRESSION 

o Simple to calculate natural frequencies, 

DISPLACEMENTS AND STRESSES DUE TO 
HARMONIC AND RANDOM FORCED VIBRATION 
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PAAM INPUT INCLUDES 

STATISTICAL DESCRIPTIONS OF RANDOM VARIABLES 


Input Variables for Lox Post Thick Cylinder Model 


Variable 

Distribution 

Mean 

COV 

r, 

truncated Normal 
(^0.003) 

0 . 094 in. 

1 .06% 


truncated Normal 
(-0.002, 40.01) 

0.110 in. 

U. 5 r j% 

r 

Normal 

3.40E+7 psi 

2% 

V 

Normal 

0 . 3594 

2% 

a 

Normal 

5.65E-6 /R 

5% 

Pi 

Lognormal 

3077 psi 

4% 

IK 

Lognormal 

3232 psi 

4% 

T, 

Lognormal 

19A R 

1.55% | 

__ 

Lognormal 

lAVi R 

1.5b% I 


o User supplies mean value, 

DISTRIBUTION TYPE, AND COV 
FOR EACH RANDOM VARIABLE 

o Many distribution types 

POSSIBLE (INCLUDING TRUNCATED) 

o Set C0V=0 for 

DETERMINISTIC VARIABLES 

0 If STATISTICAL INFORMATION 
UNKNOWN, MAKE ESTIMATES 
FROM DEFAULT SPECIFICATIONS 


PAAM OUTPUT INCLUDES CDF OF RESPONSE FUNCTION 


O FPI CALCULATES CUMULATIVE 
DISTRIBUTION FUNCTION (CDF) 
OF RESPONSE VARIABLE 


0 CDF BASED ON EQUIVALENT 
NORMAL DISTRIBUTION 
OF RESPONSE VARIABLE 


o Excellent agreement between 
FPI and Monte Carlo 
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o FPI much faster than Monte Carlo 


URiQINAL PAGE IS 
OF POOR QUALITY 
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PAAM CALCULATES SENSITIVITY FACTORS OF RANDOM VARIABLES 


LOX POST (THICK CYLINDER) 


*3 Stmdird Doutien* 



JUXDOil VARUBLI 


o Shows relative influence of 

RANDOM VARIABLES ON THE TOTAL 
UNCERTAINTY OF THE RESPONSE 

o Function of both 

- MATHEMATICS OF THE EQUATION 

- DISTRIBUTION TYPE AND COV 

OF THE RANDOM VARIABLES 

o Helps the designer/analyst to 

DETERMINE REQUIREMENTS FOR 
FURTHER DESIGN/ANALYSIS 

- MORE ACCURATE INPUT STATISTICS 

- MORE EXACT ANALYSIS MODEL 

- TIGHTER TOLERANCES IN DESIGN 

- BEST OPPORTUNITIES FOR 

DESIGN IMPROVEMENTS 


SUMMARY 


O PAAM MAKES A RAPID, INEXPENSIVE ESTIMATE OF THE 
PROBABILISTIC STRUCTURAL RESPONSE OF A SYSTEM 

o Requirements: 

- SIMPLIFIED REPRESENTATION OF COMPONENT GEOMETRY AND LOADS 

- CLOSED- FORM SOLUTIONS DERIVED OR ADAPTED FROM THE LITERATURE 

- STATISTICAL DISTRIBUTIONS OF INPUT VARIABLES 

o Probabilistic calculation performed with FPI algorithm 
o Output includes: 

- CDF OF RESPONSE VARIABLE 

- SENSITIVITY FACTORS FOR RANDOM VARIABLES 

0 PAAM PARTICULARLY USEFUL IN EARLY STAGES OF DESIGN/ANALYSIS 
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